Advanced resistivity model for arbitrary magnetization orientation applied to a series 
of compressive- to tensile-strained (Ga,Mn)As layers 
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The longitudinal and transverse resistivities of differently strained (Ga,Mn)As layers are theo- 
retically and experimentally studied as a function of the magnetization orientation. The strain in 
the series of (Ga,Mn)As layers is gradually varied from compressive to tensile using (ln,Ga)As tem- 
plates with different In concentrations. Analytical expressions for the resistivities are derived from 
a series expansion of the resistivity tensor with respect to the direction cosines of the magnetization. 
In order to quantitatively model the experimental data, terms up to the fourth order have to be 
included. The expressions derived are generally valid for any single-crystalline cubic and tetrag- 
onal ferromagnet and apply to arbitrary surface orientations and current directions. The model 
phenomenologically incorporates the longitudinal and transverse anisotropic magnetoresistance as 
well as the anomalous Hall effect. The resistivity parameters obtained from a comparison between 
experiment and theory are found to systematically vary with the strain in the layer. 
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I. INTRODUCTION 



The implementation of ferromagnetism in III-V semi- 
conductors by incorporating high concentrations of mag- 
netic elements into the group-Ill sublattice has opened 
the prospect of extending conventional semiconductor 
technology to magnetic applications i^i^i^. A prominent ex- 
ample, which is presently intensely studied, is the di- 
luted ferromagnetic semiconductor (Ga,Mn)As. Even 
though the Curie temperatures reported so far are well 
below room temperature, it represents a potential can- 
didate or at least an ideal test system for spintronic 
applications due to its compatibility with the stan- 
dard semiconductor GaAs. During the last decade, 
considerable progress has been made in understanding 
the basic structural, electronic, and magnetic proper- 
ties of (Ga,Mn)Asi^ Longitudinal anisotropic magne- 
toresistance (AMR)ii^i^il>2. and transverse AMR, often 
called planar 



feet (AHE) 



Hall effect (PHE)^ anomalous Hall ef- 



10.11.12 



and magnetic anisotropy (MA)i^i 



13.14 



have been established as characteristic features, making 
(Ga,Mn)As potentially suitable for field-sensitive devices 
and non-volatile memories i^ii^i^i Great effort has been 
made to understand the microscopic mechanisms behind 
the observed magnetic phenomena and to obtain, theo- 
retically and experimentally, values for the correspond- 
ing physical parameters. In particular, the storage and 
processing of information by manipulating the magne- 
tization as well as the readout via electrical signals de- 
mand a precise knowledge of the parameters controlling 
the MA and the AMR. It has been shown that both 
the MA and the AMR are affected by temperature, hole 
density, and strain For example, chang- 
ing the strain from compressive to tensile, the surface 
normal which usually represents a magnetic hard axis in 
compressively strained (Ga,Mn)As layers turns into an 
easy axis in tensile-strained layersJ^i^ 



As demonstrated in Ref. 25, angle-dependent mag- 
netotransport measurements, performed at different 
strengths of an external magnetic field, are a gen- 
uine alternative to ferromagnetic resonance (FMR) 
spectroscop y^^'^^ for probing the MA in (Ga,Mn) As. The 
application of this method, however, requires analytical 
expressions for the longitudinal and transverse resistivi- 
ties piong and ptrans, respectively, correctly describing the 
AMR and the AHE. Based on symmetry considerations, 
such expressions can be derived in a phenomenological 
way by writing the resistivity tensor as a series expan- 
sion with respect to the direction cosines of the magneti- 
zation. In Ref. [25I the angle-dependent magnetotransport 
data of compressively strained (Ga,Mn)As layers, grown 
on (001)- and (113)A-oriented GaAs substrates, could be 
well simulated considering only terms up to the second 
order. 

In the present work, we systematically study the in- 
fluence of vertical strain on the AMR and the AHE, 
investigating a series of compressive- to tensile-strained 
(Ga,Mn)As layers, grown on (In,Ga)As templates with 
different In contents. An advanced macroscopic model 
is presented which phenomenologically describes the de- 
pendence of piong and ptrans On the magnetization orien- 
tation for cubic ferromagnets with tetragonal distortion 
along [GDI]. The analytic expressions for piong and ptrans 
are first discussed for a variety of configurations with in- 
plane and out-of-plane magnetization and are then used 
to analyze the angle-dependent resisivities recorded from 
the (Ga,Mn)As layers under study. In contrast to Ref. [25l. 
distinct features in the longitudinal out-of-plane AMR 
occured which can only be described by taking into ac- 
count terms up to the fourth order in the magnetiza- 
tion components. Finally, the resistivity parameters, de- 
termined by fitting the calculated curves to the experi- 
mental data, are discussed as a function of the vertical 
strain in the layer. The model presented in this work ap- 
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plies not only to (Ga,Mn)As but most generally to single- 
crystalline cubic and tetragonal ferromagnets. 



II. EXPERIMENTAL DETAILS 

A series of differently strained (Ga,Mn)As layers with 
constant thickness of 180 nm and Mn concentration 
of ^ 5% was grown by low-temperature molecular-beam 
epitaxy on (In,Ga)As templates in the following way: 
After thermal deoxidation, a 30 nm thick GaAs buffer 
layer was grown at a substrate temperature of Tg « 
580 °C on semi-insulating GaAs(OOl). Then the growth 
was interrupted, Ts was lowered to '^430°C, and a 
graded (In,Ga)As layer with a thickness between /xm 
and 5 /xm was deposited following the method described 
in Ref. [27l In order to minimize the number of thread- 
ing dislocations and to end up with different lateral lat- 
tice constants in the (In,Ga)As templates, the In con- 
tent was continuously increased in each template from 
2% up to a maximum value of 13%. Prior to the epi- 
taxy of (Ga,Mn)As, the growth was again interrupted 
and Ts was lowered to ~250°C. High-resolution x-ray 
diffraction (HRXRD) reciprocal space mapping (RSM) 
of the (224) reflex was used to determine the vertical 
strain Ezz — {a± — arci)/arci of the (Ga,Mn)As layers, 
where the relaxed lattice constants Orci were derived 
from the lateral and vertical lattice constants a\\ and 
a^, respectively, applying Hooke's law. The values of 
Ezz were found to gradually vary from +0.24% for the 
compressive-strained sample without (In,Ga)As template 
to —0.46% for the tensile-strained sample with 13% In. 
Moreover, RSM showed that the (In,Ga)As layers were 
almost completely relaxed whereas the (Ga,Mn)As layers 
were fully strained. Further details of the growth proce- 
dure and the RSM method will be presented elsewhere. 

For the magnetotransport studies two types of Hall 
bars with current directions along [100] and [110] were 
prepared on several pieces of the cleaved samples. The 
width of the Hall bars is 0.3 mm and the longitudi- 
nal voltage probes are separated by 1 mm. High-field 
magnetotransport measurements (up to 14.5 T) at 4.2 K 
yielded hole densities for the as-grown samples between 
3 X lO^'^cm"^ and 4 x 10^°cm~^. Least squares fits 
were performed to separate the contributions of the nor- 
mal and anomalous Hall effect. Curie temperatures be- 
tween 61 K and 83 K were estimated from the peak po- 
sitions of the temperature-dependent sheet resistivities 
at 10 mT. For the angle-dependent magnetotransport 
measurements, carried out at 4.2 K, the Hall bars were 
mounted on the sample holder of a liquid- He-bath cryo- 
stat which was positioned between the poles of an elec- 
tromagnet system providing a maximum field strength 
of 0.68 T. The sample holder has two perpendicular axes 
of rotation, allowing for any orientation of the Hall bars 
with respect to the applied magnetic field H . 



III. THEORETICAL MODEL 

The macroscopic theoretical model presented in this 
paper is based on the assumption that the sample 
area probed by magnetotransport can be approximately 
treated as a single homogeneous ferromagnetic domain. 
It provides analytical expressions for the electrical resis- 
tivities as a function of the magnetization orientation. 
Although the single-domain picture is known to usually 
break down in situations where the magnetic system is 
undergoing a magnetization reversal process, it has been 
successfully applied to the description of a variety of 
magnetization-related phenomena in (Ga,Mn)As, partic- 
ularly at sufficiently high external magnetic fields. 

Given a single ferromagnetic domain, the macroscopic 
magnetization is described by the vector M — Mm 
where M denotes its magnitude and the unit vector m 
with the components m^;, TOj,, and its orientation. 
Throughout this paper, all vector components labeled by 
X, y, and z refer to the cubic or tetragonal coordinate 
system associated with the [100], [010], and [001] crystal 
directions, respectively. 



A. Longitudinal and transverse resistivities 

We consider a standard configuration for magneto- 
transport measurements where the longitudinal voltage 
is probed along and the transverse voltage across the di- 
rection ji of a homogeneous current with density J = Jj. 
Accordingly, we introduce a right-handed coordinate sys- 
tem with unit vectors j, and n = j xt, so that t defines 
the transverse direction and n typically the surface nor- 
mal. The measured voltages arise from the components 
-E-iong = j ■ E and i?trans — t ■ E the electric field E. 
Starting from Ohm's law E ~ p ■ J with the resistivity 
tensor p, the longitudinal resistivity piong (sheet resistiv- 
ity) and transverse resistivity ptrans (Hall resistivity) can 
be written as 



long 



St, 



3-P-3 



tpj- (1) 



1. Resistivity tensor 

Following the ansatz of Birss,^® the dependence of piong 
and ptrans On the magnetization orientation m is derived 
in a phenomenological approach by writing the compo- 
nents pij of the resistivity tensor p as series expansions 
with respect to TOx, ruy^ and niz. In Ref. [2^ we have pre- 
sented analytical expressions for piong and ptrans includ- 
ing terms up to the second order. Now, forced by the ex- 
perimental results, we substantially extend this model by 
taking into account terms up to the fourth order. Using 
the Einstein summation convention, the series expansions 
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read as 



'Hj + dkij-rrik + akii-jnikmi + akimijmkmimm - 



(2) 

where the components of the galvanomagnetic tensors 
o-ijj o-kiji ■■■ appear as expansion coefhcients. Neumann's 
principle requires that a tensor representing a macro- 
scopic physical property of a crystal must be invariant 
under all symmetry operations S of the corresponding 
point group. The mathematical formulation of this re- 
quirement leads to conditional equations for the expan- 
sion coefficients reading as 



ttij 
a-kiij 

^klmnij 



\S\SkqSioS jpdqop , 
\S\ ^kqSlrSioS jpdqrop i 
\^\ SkqSlrSjYitSioS jpClqrtop •> 
\S\ SkqSlrSmtSnuSioS 



jpCLqrtuop ; 



(3) 



where Stj denotes the ij component and \S\ the deter- 
minant of the symmetry matrix S. The determinant ap- 
pears in the formula of Eqs. ^ since M is an axial vec- 
tor, whereas E and J are polar vectors. In order to derive 
a complete set of conditional equations for the expansion 



coefficients, it is sufficient to apply Eqs. ((S]) to a small 
set of generating symmetry matrices. In the cases of cu- 
bic symmetry Td and tetragonal symmetry D2d this set 
consists of only two matrices, namely Ss,Sg and 52,5*8, 
respectively. The matrices are given bjj2^ 



52 = 



-1 








-1 




1 


• 


, 58= 1 







• 









[0 








So = 




(4) 



Insertion of the generating matrices into Eqs. ([3]) reveals 
that most of the expansion coefficients are equal in pairs 
or vanish. The resulting resistivity tensor for tetragonal 
symmetry can be separated into two terms 



Ptctragonal — Pcubic 



Ap 



(5) 



where pcubic is the resistivity tensor for cubic symmetry 
and Ap a difference term which vanishes in the case of 
perfect cubic symmetry. Written in ascending powers of 
nix, my, and m^, the two parts of ptctragonai are given by 



Pcuhu 




m^myTn"^ 
rnxXn^mz m'^mymz 




rnxmyiriz 
mlmymz 



m^my 






(6) 



A-p 




e4ml 
esm^myml 





02171x771;, 
,2 

A 



C27nxmz 
CzTnl C27ny7nz 
C27nymz Ciml 




esmxTnyTTi, 
e2TO^m^ -I- 6^7711 



ee77ix7n^7nz 

2 

eQmx77ly77lz 



■ eiTTlxTTlz 

■ erTUyTnl 

4 

- eim^ 





-d277lxm1 



—d27ny7nl 
d27nx7n1 




(7) 



Deriving the above equations, we repeatedly made use of 
the trivial identity 



|mp = 7nl+my + 7nl = l. 



(8) 



The expansion parameters A, E^, and a, ey are 
linear combinations of the non- vanishing expansion coef- 



ficients, reading as 

A 
B 
Ci 

C2 

D 



El 
E2 
E3 



flu + aii22 + 0111122 , 

0-123 + 3ai2223 : 

aim — 11122 + 6aii22ii — 2aiiii22 , 
2ai2i2 + 4aiii2i2 , 

aill23 — 3ai2223 , 

aiiiiii 6aii22ii -I- aiiii22 , 
60112233 — 2aiiii22 , 

12aii2323 ~ 4aiii2i2 • 



(9) 
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a = 033 ^ ail + flll33 ~ ail22 + flllll33 ^ 0,111122 , 

b = 0312 — ai23 + 3aii3i2 — 8012223 j 
Ci = 03333 ^ aim + 01122 — aii33 + 6aii3333 — 60112211 
+ 2aiiii22 — 2aiiii33 1 

Cl — 2a2323 — "iO-YlVl + 4aiii3i3 — 4aiii212 , 



C3 


— 03311 — 11122 + 6aii33ii — 6aii2211 , 


rfl 


— 3ai2223 — 


3ail312 + 033312 ~ 011123 7 


rf2 


= 3ai3323 — 


3ai2223 , 


ei 


= 0333333 — 


oiiiiii + 6aii22ii — 60113333 




+ aiiii33 


— 0111122 , 


62 


= 6aii3322 - 


- 60112233 + 60112211 — 60113311 


63 


— 12ai233i2 


— 12oii2323 I 


64 


— 0333311 ~ 


0111122 + 60112211 — 60113311 , 


65 


= 2aiiii22 - 


- 2oiiii33 1 


66 


= 4aiii2i2 - 


- 4oiii3i3 , 


67 


= 4ai333i3 - 


- 4oiii3i3 . 



It should be noted that the present notation is partiaUy 
different from that used in Ref. 25. For cubic ferromag- 
nets with a small tetragonal distortion along [001], as in 
the case of the (Ga,Mn)As layers under investigation, the 
expansion parameters are expected to linearly vary with 
Ezz- The parameters A and a, just as all other expansion 
parameters, then read as 

Arel+£zzA! , a^Ezzo! . (11) 

While o vanishes for zero strain, A becomes identical to 
Arex describing the resistivity of the relaxed cubic crystal. 

Once the resistivity tensor p is known, Eqs. ([1]) allow 
to calculate piong and ptrans for any current direction j 

I 



and any orientation t of the transverse voltage probe rel- 
ative to the crystal axes. For a concise presentation of 
the results, it is convenient to replace the components 
Tfix, ruy, and of m referring to the cubic or tetrag- 
onal crystal axes with mj, rrit, and to„ referring to the 
more experiment-related coordinate system defined by j, 
t, and n according to 

m, = j^rrij + Umt + niTTin , {i = x,y,z). (12) 

Consequently, Eq. ([8]) has to be rewritten as 

Imp = + + to2 = 1 . (13) 

For the rest of the paper we exclusively focus on the most 
common case of a current flow parallel to the surface of 
a (OOl)-oriented sample and examine in detail the two 
situations where the current direction is along [100] and 
[110]. 

2. piong an-d ptrans fov Current in the (001) plane 



To cover first the general case of a current flowing along 
an arbitrary direction within the (001) plane, we write j, 
t, and n as 




where a denotes the angle between the current direction 
and the [100] crystal axis. Using Eqs. H]), and 
the resistivities piong and ptrans can be written as poly- 
nomials of fourth order in the variables rrij , mt , and m„ : 



Plo 



Ptr 



A+[C,-C, + \Ey-^ 



1 — cos 4q; 



Cl - (Cl - C2 - El 



El cos 4a mj 



-E2 + 64 - 



-El + E2 + E2 



C3 + E2- [Cl -C2+E1+E2+E3 
1 — cos 



1 — cos Aa 



E2 + { El — E2 — E:- 



1 — cos 



— (Cl — C2) sin4Q! rrijmt — Ei sin 4a m'jmt + — ^£^2 + E3 ) sin 4a mjUitrri^ , 



2 

-BtJln 



C2 + [Cl ^ C2 + El) 



1 — cos 4a 



rrijmt — Dmz 



E3- E1+E2 + E, 



1 — cos 4a 



mjmtm^ 



i (Cl -C2 + El) sin 4a {ml - mJ) + i [Ei +E2 + -E3) sin 4a (to| - mj) mri^ 



(15) 



(16) 
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For simplicity, we have introduced the abbreviations 
B = B + b , D = D + di , 



E2 



62 



E-l 



63 



(17) 



Equations ([TS]) and drasticaUy simphfy for cur- 
rent directions along [100] and [110], referring to a = 0° 
and a = 45°, respectively. In both cases the correspond- 
ing expressions for the resistivities take the same form, 
namely 



Plong 



Ptrans = PQrrin + pTrnjiTit + psm^ 

2 



(18) 



(19) 



-I- pgrrijiTitni^ 

The parameter po in Eq. ^TE\i may be regarded as a ref- 
erence for the angular dependence of piong- It represents 
the longitudinal resistivity for parallel alignment between 
m and t where mj = nin = 0. The reference direction 
t can be easily changed to j or n by substituting rrij or 
nin, respectively, with the help of Eq. (fT3|l . All resistiv- 
ity parameters pi {i = 0, 9) are linear combinations 
of the expansion parameters A, B, Ci, and depend, 
according to Eqs. (fTSjl and (fT6|) . on the angle a, i.e., on 
the direction j of the current with respect to the crystal 
axes. 

For j II [100] and t \\ [010] the resistivity parameters pi 
are given by 



(20) 



For unstrained layers with perfect cubic crystal symme- 
try the expansion parameters represented by small letters 
vanish and we obtain the relation P4 = Ps = —p2- Ac- 
cordingly, Eq. (fTSjl reduces to 



Po 


= A, 


P5 


- -E2 


Pi 


= Ci , 


P6 


= -B, 


P2 


= + C3 , 


Pi 


= C2, 


P3 


= ^L' 


Ps 




Pi 


= -£'2 -f- 64 , 


P9 


— E3 . 



Pi 



ong 



Po 



Pim] 



P2m1mi 



(21) 



For j II [110] andt 
are given by 



Po 


= A 




2[^^ 


Pi 


= C2 


+ 


El, 




1 






P2 


2 




Ci + 


P3 


= ~Ei 






1 


"1 




Pi 


2 


2 


El - 


P5 


- -El 


+ E3 


P6 


= -B, 




P7 


= Ci 


+ 


El , 


P8 


= -D, 




Po 


= -El 


— E2 



[110] the resistivity parameters pi 



1 



C2] + ::Ei, 



C2 + 2c3 — £/i + E2 — E^ 



E3 + 264 



(22) 



Here, perfect cubic crystal symmetry yields the con- 
straints 2p2 = Pi + 3p3 — P5 — P7 ~ P9 and 4p4 — 
-5p3 + 2p5 + 2pg. 

The parameters pi {i — 0, 9) may be thought of as 
the components of ten-dimensional vectors p which for 
j II [110] and j II [100] are related to each other by the 
linear transformation 



P[iio] = T ■ P[ioo] • 
The matrix T is given by 



(23) 



T = 



1 
4 

. 1 

1 -k 
. -1 

1 
4 



1 



1 



1 



(24) 



■ .-1.1. 

and satisfies the equation T^^ ~ T. 

3. Polar plots of piong and ptrans 

In order to graphically illustrate the dependence of 
Piong and Ptrans On the magnetization orientation m, it is 
instructive to set either m„, ruj, or rut in Eqs. and 
p9)) equal to zero and to display the resulting analytical 
expressions for piong and ptrans in polar plots. 

For in-plane magnetization the identity 77i„ = holds 
and Eqs. HH]) and simplify to 



Pim] 



Piong — Po 
Ptrans = PrmjUlt 



Psnij 



(25) 
(26) 



In Fig. [T] the angular dependence of piong — Po is depicted 
by the solid line, clearly reflecting the longitudinal in- 
plane AMR, i.e., the variation of piong with the magneti- 
zation orientation. The dashed line illustrates the angu- 
lar dependence of ptrans and reflects the transverse AMR. 
The dotted circle is the separating line between negative 
and positive values of piong — Po and ptrans- 

An orientation of the magnetization perpendicular to 
the current direction is equivalent to the condition mj — 
0. In this case Eqs. and ([TO]) reduce to 



Piong — Po 
Ptrans — P6"^n 



P2"4 



P4TO* 



psml 



(27) 
(28) 



The corresponding polar plots of piong — Po and ptrans 
in Fig. [2] reflect the longitudinal out-of-plane AMR for 
m L j and the AHE, respectively. In the special case of 
cubic crystal symmetry and j \\ [100], Eq. (P7)) further 
reduces to 



Plo 



Po 



P2mtml , 



(29) 
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FIG. 1: Schematic polar plot in arbitrary units showing the 
dependence of piong — po (sohd line) and ptrans (dashed line) 
on m for in-plane magnetization perpendicular to n. The 
dotted circle indicates the zero line which separates negative 
and positive values. 

according to Eq. (|2T|) with nij = 0. As demonstrated by 
the dash-dotted line in Fig. [21 piong now exhibits fourfold 
symmetry. 




out-of-plane AMR) 

FIG. 2: Dependence of piong— po (solid line) and ptrans (dashed 
line) on m for out-of-plane magnetization perpendicular to 
j. The dash-dotted line illustrates the special case of cubic 
crystal symmetry and current along [100]. 

The case m Lt, or equivalently rrit = 0, is illustrated 
in Fig. [3l where the solid line again describes the longi- 
tudinal out-of-plane AMR and the dashed line the AHE. 
The appropriate equations for piong and ptians read as 

Piong = Po + Pirn^j + P2ml + psm'^ + p4m'^ 

+ psm^^ml , (30) 

Ptrans = Psmn + Psm^ . (31) 

The expression for ptrans in Eqs. ((^ and ([31]) prop- 
erly describes the AHE in a phenomenological way, but 




out-of-plane AMR) 

FIG. 3: Schematic polar plot in arbitrary units showing the 
dependence of piong — po (solid line) and ptrans (dashed line) 
on m for out-of-plane magnetization perpendicular to t. 

it does not include the ordinary Hall effect. For mag- 
netic field strengths iiqH < 1 T and hole concentrations 
p > 10^^ cm~^ as in our experiments, however, the contri- 
bution of the ordinary Hall effect noH/ep (e denotes the 
elementary charge) to ptrans is smaller than 6 x 10~^ cm 
and thus about two orders of magnitude smaller than the 
measured peak values of ptrans (see Section IV). Effects 
correlated with the magnitude B of the magnetic induc- 
tion B, such as the negative magnetoresistance, can be 
taken into account by considering i?-dependent resistiv- 
ity parameters. 

4- Poly crystalline materials 

In the literature, the longitudinal in-plane AMR and 
the transverse AMR in (Ga,Mn)As are often theoretically 
described by the equation a^^'^" 

Piong = p_L + (p|| - p±) COS^(/) , (32) 
Ptrans = (p|| - P_l) siu (/) COS , (33) 

where (p denotes the angle between j and m. Taking 
into account the relations cos^ (j) — "^j ^nd sin (j) cos 4> — 
rrijmt, which are only valid for an in-plane configuration, 
it becomes clear that Eq. ([32]) is a good approximation 
to Eq. (j25p only if the fourth-order term p^rri^ is neg- 
ligibly small compared to pim^, i.e., if jpaj ^ \pi\, or 
more generally ^ |Ci|,|C2|. The results presented 
in Sec. IV B reveal that these inequalities do not apply 
to any of the (Ga,Mn)As samples under study. On the 
other hand, the angular dependence of ptrans described 
by Eq. ([55)) agrees with that given by Eq. ([^ . The pref- 
actors of the cos^ (j) and sin ^ cos ^ terms, however, are in 
general different and equal only in the special case where 
Pi = pj, or equivalently Ci = C2. As already pointed 
out in Ref. i25|, Eqs. ([5^ and (|33p only apply to isotropic 
materials such as polycrystals where the resistivities do 
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not depend on the direction j of the current relative to 
the crystal axes. In fact, Eqs. (|32p and ([55]) result from 
Eqs. (HJ) and ^ by averaging over all possible spatial ori- 
entations of the coordinate system {j,t,m} with respect 
to the cubic crystal axesi^ Calculation yields 

P± = A + i (Ci - C2) + 4 i^Ei + 3E2 - E3) , 

Pl\ -P^ = I (2^^i + 3C2) + ^ il2Ei - 2E2 + iEs) . 

(34) 



the thickness of the (Ga,Mn)As sample is small, '^^ or a 
uniaxial in-plain strain is induced by a piezoelectric actu- 
ator mounted onto the sample.— 1^ The incorporation of 
the contribution Cjj cos 2^/; from Eq. ([55)1 into our model 
would yield a further term in the expression of piong in 
Eq. (fT8|) . For j \\ [100] it reads puirijmt and leads to an 
additional angular dependence of piong- For j || [110] it 
is given by — Pu/2 -I- p„m| and can be incorporated into 
the terms po and piirij. 

B. Magnetic anisotropy 



5. Uniaxial [110] in-plane anisotropy 

Recently, Rushforth et al. presented an experimental 
and theoretical study on the components of the in-plane 
AMR based on the equations^^ 



Pi 



ong 



Pa 



PtT 



Ci COS 2<j) + Cu COS 2il} + Cc cos ^ip 
Ci,c cos(4?A - 20) , (35) 
Ci sin 20 - C/,c sin(4?/' - 20) , (36) 



which were obtained by extending the model of Doring,^'^ 
introduced for cubic Ni, to systems with cubic [100] plus 
uniaxial [110] anisotropy. Here again denotes the an- 
gle between j and m, and -0 the angle between m and 
the [110] crystal direction. The coefficients C/, Cu, Cc, 
and C/.c represent a non-crystalline, a uniaxial, a cu- 
bic, and a crossed non-crystalline/crystalline contribu- 
tion, respectively, pav is the average value of piong as M 
is rotated through 360°. Starting from Eqs. (|15p and pB]) 
with TO„ = 0, it can be shown that the coefficients C/, 
Cc, Ci^c, and pav are related to the parameters A, Ci, 
C2, and El, introduced in Sec. Ill A, by the equations 



Ci 
Cc 
Ci^c 

Pav 



1 



^Pav 
1 



S>Pav 
1 



(Ci +C2 + El) , 



El , 



(Ci — C2 + El 



A+lci + '-Ei 



(37) 



In our model no uniaxial [110] anisotropy of the resis- 
tivity has been taken into account, i.e., Cu — 0, since 
no significant influence of such a contribution has been 
found in analyzing the angle-dependent magnetotrans- 
port data. This is consistent with the fact that the 
strength of the magnetic uniaxial in-plane anisotropy (for 
definition see Sec. Ill B), inferred from the experimental 
data, turned out to be negfigibly small in the (Ga,Mn)As 
samples under study. A uniaxial in-plane contribution 
may become important, however, when the experiments 
are performed at temperatures much higher than 4.2 K,^^ 



In the preceeding section analytical expressions for 
Piong and ptrans as a function of the magnetization ori- 
entation m have been derived. The direction of m, in 
turn, is determined by the MA of the ferromagnetic ma- 
terial and by the strength and orientation of an external 
magnetic field H. 

Magnetic anisotropy stands for the dependence of the 
free energy density _F of a magnetic system on the orien- 
tation of M . In addition to supposing a simple single- 
domain model, we assume that the magnitude M of the 
magnetization is nearly constant under the given experi- 
mental conditions. Instead of F we therefore consider the 
normalized quantity Fm = F/M , allowing for a more con- 
cise description of the MA. For a (Ga,Mn)As film with 
tetragonal distortion along [001], the anisotropic part of 
Fm can be written a o-^^i^^ 



FM{rn) = (m'^ + my) + Bi^m'l + B2i_m\ 



PoM 



+ Biioim^ - TTiyf 



(38) 



The first three terms are intrinsic contributions arising 
from spin-orbit coupling in the valence band. The fourth 
and fifth terms are extrinsic contributions describing 
the demagnetization energy of an infinite plane (shape 
anisotropy) and a uniaxial in-plane contribution, respec- 
tively, whose origin is still under discussion.^ '^'^i'^^ The 
two ml terms cannot be distinguished in our experiments 
and are therefore lumped into a single term BoQim^^. 

In the presence of an external magnetic field H one 
has to additionally take into account the Zeeman energy, 
and the total energy density is finally given by the (nor- 
malized) free enthalpy density 



Gji/(m) = FM[m) 



Mo 



H 



(39) 



Figure m shows a graphical illustration of the various con- 
tributions to Gj\/. Given an arbitrary magnitude and 
orientation of the magnetic field H , the direction of the 
magnetization M is determined by the minimum of Gm 
with respect to the components of m. 



IV. RESULTS AND DISCUSSION 

For sufficiently high magnetic fields the contribution of 
the free energy Fm to the free enthalpy Cm in Eq. 
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Bc4||(m,+ my)+B^,_Lm^ B„„, By 




Ai(,H m 



FIG. 4: Graphical illustration of the different contributions 
to the normalized density of the free enthalpy Gm, plotted as 
three-dimensional functions of the magnetization orientation. 
The equilibrium position of M is determined by the minimum 
of Gm ■ 



becomes much smaller than the contribution of the Zee- 
man energy and can, as a good approximation, be ne- 
glected. In this case the magnetization M aligns with H 
and the dependence of piong and ptrans on the magnetiza- 
tion orientation can be simply probed by systematically 
varying the direction of H. The values of the resistiv- 
ity parameters pi are then derived by fitting Eqs. (jlSp 
and to the measured data with the vector compo- 
nents of m replaced by those of h. When the magnetic 
field is gradually lowered, however, the Zeeman term in 
Eq. p9|l decreases, the relative contribution of Fm to 
Gm increases, and the orientation of M more and more 
deviates from the direction of H towards one of the easy 
axes determined by the minima of Fm- In other words, 
the measured resistivities are increasingly influenced by 
the MA. In Ref. [2^ we have shown that this effect can 
be utilized to probe the MA in (Ga,Mn)As by means of 
angle-dependent magnetotransport measurements. 

In our experimental setup the field strength of the elec- 
tromagnet is limited to 0.68 T. In most cases this value 
suffices to align M almost perfectly along H. In situ- 
ations, however, where H approaches a hard magnetic 
axis, this might no longer apply and the influence of the 
MA has to be taken into account when deriving resistiv- 
ity parameters from angle-dependent resistivity curves. 
Therefore, we routinely determined both resistivity and 
anisotropy parameters for all samples under study. The 
corresponding procedure is exemplified in Sec. IV A by 
means of an almost unstrained (Ga,Mn)As layer with Szz 
= —0.04%. In Sec. IV B the strain dependence of the re- 
sistivity parameters pi will be discussed. An analogous 
study on the anisotropy parameters Bi goes beyond the 
scope of this paper and will be presented elsewhere. A 
brief survey, however, is given in Ref. 



A. Determination of the resistivity parameters 

The longitudinal and transverse resistivities of the 
(Ga,Mn)As layers were measured for both j \\ [100] and 
j II [110] as a function of the magnetic field orientation at 
fixed field strengths of hqH = 0.11, 0.26, and 0.65 T. At 
each field strength, H was rotated within three diff'erent 
crystallographic planes perpendicular to n, j, and t, re- 
spectively. The corresponding configurations, labeled I, 
II, and III, are shown in Fig. \5\ In the case of M \\ H 
they are identical to the configurations used in Figs. [TH3] 
to illustrate the theoretical angular dependences of piong 
and ptrans by means of schematic polar plots. 




FIG. 5: The angular dependence of the resistivities was 
probed by rotating an external magnetic field H within three 
different planes (a) perpendicular to n, (b) perpendicular to 
j, and (c) perpendicular to t. The corresponding configura- 
tions are referred to as I, II, and III. 

Prior to each angular scan, the magnetization M was 
put into a clearly defined initial state by raising the field 
to its maximum value of 0.68 T where M is supposed 
to nearly saturate and to align with the external field. 
The field was then lowered to one of the above men- 
tioned magnitudes and the scan was started. Figures [6] 
and [7] show as an example the angular dependences of 
Piong and ptrans for the almost unstrained sample with 
£zz — —0.04%, measured for j along [100] and [110], re- 
spectively. The experimental data are depicted by red 
circles. In the first case, H was rotated in the (001), 
(100), and (010) plane, in the second case the planes of 
rotation were the (001), (110), and (110) plane, corre- 
sponding to the configurations I, II, and III, respectively. 
At 0.65 T the Zeeman energy dominates the free enthalpy. 
Consequently, M nearly aligns with H and continuously 
follows its motion. In fact, the curves of piong and /Otians 
at 0.65 T are smooth and largely reflect the angular de- 
pendences of the resistivities described by Eqs. (fT5|) and 
(fro| with m replaced by h. With decreasing magnetic 
field, the influence of the MA increases and the orienta- 
tion of M deviates more and more from the fleld direc- 
tion. Accordingly, jumps and kinks occur in the curves 
at 0.26 and 0.11 T, arising from sudden movements of M 
caused by discontinuous displacements of the minimum 
of Gm. 

Values for the resistivity parameters pi and the 
anisotropy parameters Bi were determined by an iter- 
ative fit procedure. Starting with an initial guess for the 
anisotropy parameters, the resistivity parameters were 
obtained by fitting Eqs. (fT5|) and ITOl) to the experimen- 
tal data recorded at 0.65 T. Then the anisotropy parame- 
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FIG. 6: (Color online) Resistivities piong and ptrans recorded from the nearly unstrained (Ga,Mn)As layer with Szz = —0.04% 
at 4.2 K and j \\ [100] (red circles). The measurements were carried out at fixed field strengths of fj.oH = 0.11, 0.26, and 0.65 T 
with H rotated in (a) the (001), (b) the (100), and (c) the (010) plane, corresponding to the configurations I, II, and III, 
respectively. The black solid lines are fits to the experimental data using Eqs. (|18|l and p9|l . and one single set of resistivity 
and anisotropy parameters. The dashed line in (b) at 0.65 T represents an attempt to fit the measured curve considering only 
terms up to the second order. 



TABLE I: Values of the resistivity parameters in units of 
10~*fl cm determined for the nearly unstrained (Ga,Mn)As 
layer with e^^ = -0.04%. 



p^ (10"" ncm) 


j II [100] 


3 II [110 


po (0.65 T) 


103.9 


103.8 


po (0.26 T) 


104.8 


104.7 


po (0.11 T) 


105.2 


105.1 


Pi 


-5.3 


-2.3 


P2 


-2.2 


-1.9 


P3 


1.7 


-1.7 


Pi 


2.2 


2.2 


P5 


2.2 


-0.4 


P6 


9.8 


9.8 


P7 


-4.0 


-3.6 


PS 


0.0 


0.0 


P9 


1.3 


0.5 



ters were modified for an optimal agreement at 0.26 and 
0.11 T, and the whole procedure was repeated until no 
further improvement of the fit could be achieved. The 
unit vector m, whose components enter Eqs. ([TSl) and 
P^ . was calculated for any given magnetic field H by 
numerically minimizing Gm with respect to the direction 
of m. With the exception of po the resistivity parame- 



ters were assumed to be field independent, which turned 
out to be a good approximation within the accuracy of 
the fit. pq was found to decrease with increasing mag- 
netic field, reflecting the negative-magnetoresistance be- 
havior of piong- The values of the resistivity parameters 
for j II [100] and j \\ [110], obtained from the fits, are 
listed in Table HI The two sets of parameters are related 
to each other according to Eqs. l|^ and ([M)) . in agree- 
ment with the theoretical model presented in Sec. III. 
The parameter pg is not immediately accessible by mea- 
surements performed in the configurations I, II, and III, 
since the corresponding term in Eq. (|19p vanishes in all 
three cases. It can be determined, however, either di- 
rectly by orienting H in a way that rrij , mt , and m„ are 
all different from zero, or indirectly by using the relation 
p'g = p^ — ps, where the primed and unprimed parameters 
correspond to j \\ [110] and j \\ [100], respectively, and 
vice versa. For the anisotropy parameters we obtained 
the values B^^ = i?4_L — — 35mT, i?ooi — 35 mT, and 
BiiQ — — 5mT. The theoretical curves calculated with 
these parameters are drawn as solid lines in Figs. [5] and 

m 

For /io-ff — 0.26 and 0.65 T, an excellent agreement 
between the measured and simulated curves is achieved, 
while for (J-qH = 0.11 T significant differences emerge. 
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FIG. 7: (Color online) Continuation from Fig.|6l The current direction is now along [110] and the magnetic field H was rotated 
in (a) the (001), (b) the (110), and (c) the (llO) plane, corresponding to the configurations I, II, and III, respectively. 



We interpret these differences as clear evidence for the 
gradual breakdown of the single-domain model with con- 
stant magnetization magnitude M at low magnetic fields. 
The experimental curves recorded in the configurations 
II and III at 0.11 T are much smoother than the calcu- 
lated ones, probably refiecting the formation of a multi- 
tude of differently oriented ferromagnetic domains. This 
interpretation is supported by a number of investiga- 
tions visualizing the domain structure of (Ga,Mn)As lay- 
ers at low magnetic fields^'^'^S'^i^ as well as by fur- 
ther magnetotransport measurements on the (Ga,Mn)As 
samples under study, not presented in this paper. Re- 
markably, the theoretical curves in Figs. [6ja) and[7ja), 
calculated for configuration I where H is rotated in the 
(001) layer plane, almost perfectly describe the measured 
curves even at 0.11 T. This may be explained by the fact 
that due to shape anisotropy the [001] axis perpendicu- 
lar to the surface is a hard axis even in the unstrained 
layer and that, as a consequence, ferromagnetic domains 
with in-plane magnetization remain more stable at low 
magnetic fields than domains with out-of-plane magne- 
tization. Accordingly, the jumps and kinks in Figs. [6]^ a) 
and[7l^a) at 0.11 T are extremely well pronounced, reflect- 
ing the occurence of a nearly perfect coherent switching 
of the whole spin system. 

Whereas the measured resistivity curves presented in 
Ref. m for two (Ga,Mn)As layers, grown on GaAs(OOl) 
and GaAs(113) substrates, could be satisfactorily well fit- 
ted by analytical expressions containing only terms up to 



the second order in the components of m, it is now com- 
pulsory to take into account higher-order terms. In order 
to illustrate the importance of the fourth-order terms for 
a correct description of the experimental traces of piong, 
simulations are depicted as dashed lines in Figs.[Sfb) and 
[TJ^b) where only terms up to the second order in were 
considered. No matter which values for p2 in Eq. (jl8|) 
were chosen, the curves completely failed to reproduce 
the measured angular dependence of piong- 



B. Strain dependence of the resistivity parameters 

It is well known that a distortion of the (Ga,Mn)As 
crystal lattice leads to a significant change of the MA. 
Microscopically, this can be explained by a strain-induced 
warping of the valence bands in addition to that caused 
by spin-orbit coupling.- The warping, however, not only 
affects the MA but also the AMR and thus the resistivity 
parameters pi. 

In order to examine the influence of the vertical 
strain Szz on piong and ptians, angle-dependent magne- 
totransport measurements analogous to those presented 
in Figs. [5] and [7] were performed on all (Ga,Mn)As sam- 
ples under study. Using the procedure described in the 
previous section, the resistivity parameters po, ps for 
j II [100] and j \\ [110] were determined independently 
of each other without taking into account their mutual 
relations theoretically predicted by Eqs. and 
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FIG. 8: Values of the normalized resistivity parameters pi/ po 
(i — 1, 8) plotted against the vertical strain Szz- The 
solid lines are smoothing splines and are drawn to guide the 
eye. In (a) the normalized values of pi, pi axe shown for 
3 II [100] and in (b) for 3 \\ [110]. The dash-dotted and dashed 
lines depict the strain dependence of pi + pa and p2 + pi, 
respectively, and are obtained by adding the individual spline 
curves. In (c) the normalized values of ps, pe + ps, and pr are 
shown for both current directions. 



Whereas the values obtained for po randomly scatter 
in the range from 5.0 x 10^'^ Q cm to 8.8 x lO^'^ cm, 
presumably due to variations in the hole concentration 
and mobility, a distinct correlation with the strain is 
found for the other resistivity parameters. This corre- 
lation is seen most clearly by considering the normalized 
quantities Pi/po instead of pi. In Fig. [5] the values of 
Pi/po {i — 1, ••■I 8) are plotted against e^z in the range 
-0.46% < < 0.22% for j \\ [100] and j \\ [110]. The 
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FIG. 9: Expansion parameters B + b, 64 plotted against 
the vertical strain Szz- The values were calculated from pi 
{i — 1, 8) using Eqs. (|20[l and (|22p . and are normalized 
to that of A. The solid lines are smoothing splines and are 
drawn to guide the eye. 



parameter pg (not shown) was determined by the indirect 
method described above and turned out to be nearly in- 
dependent of the strain with pg/po ~ 0.015 for j \\ [100] 
and p9/po«0for J || [110]. 

According to the theoretical model presented in 
Sec. Ill A, the resistivity parameters for j \\ [100] and 
j II [110] should be linearly related to each other by 
Eqs. (^5]) and yielding, for instance, the symmet- 

rical relations 



Pi P3 + P7 , 
P's = -P3 ; 

p'q = P& , 



P7 



Pi + P3 



(40) 



The primed and unprimed parameters correspond to the 
current directions [110] and [100], respectively, and vice 
versa. Inspection of Fig. [5] reveals that the measured data 
comply with the above relations supporting the validity 
of the theoretical model. 

In contrast to the resistivity parameters pi the expan- 
sion parameters B, Ci, 64, do not depend on the cur- 
rent direction and have thus a more fundamental mean- 
ing. They can be unambiguously calculated from the 
resistivity parameters using Eqs. ^U\i and (HH). If pg is 
measured directly, only one set of resistivity parameters 
is needed, otherwise the pi have to be known for both 
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current directions, j \\ [100] and j \\ [110]. In this work 
the second case apphes. Figure [3] shows the values of the 
expansion parameters normaUzed to A in dependence on 
the vertical strain £22. The smoothing spUne curves sug- 
gest that the parameters related to the cubic part Pcubic 
of the resistivity tensor (capital letters) as well as the pa- 
rameters appearing in the strain-induced difference term 
Ap (small letters) almost linearly vary with e^z in agree- 
ment with Eq. (fTTjl . 

The variation of the resistivity parameters with the 
strain manifests itself in a pronounced change of the an- 
gular dependences of piong and ptrans- This is exemplarily 
demonstrated in Figs. [TO] and [Til where the normalized 
resistivities {piong - Po) / Po and ptrans/po, respectively, at 
HqH = 0.65 T are plotted as a function of the magnetic 
field orientation and the strain for some of the config- 
urations presented in Figs. [H] and [7] The experimental 
data are again depicted by red circles and the calculated 
curves by black solid lines. In Fig. [TOTa) the longitudi- 
nal resistivity measured in configuration I is shown for 
j II [110]. Provided that M and H are appoximately 
parallel to each other at 0.65 T, the angle-dependent os- 
cillations of Piong reflect the longitudinal in-plane AMR 
illustrated in Fig. [TJ The amplitude of the oscillations is 
given by pi + p^ and increases with decreasing compres- 
sive and increasing tensile strain. It is depicted by the 
dash-dotted line in Fig. [Sfb). The specific shape of the 
oszillations is determined by the ratio of pi and p^ repre- 
senting the and m| terms, respectively. As shown in 
Fig. (IKb) , pi strongly varies with e^z whereas pa remains 
nearly constant. Figure[5Ka) reveals that for j || [100] pi, 
P3, and pi -I-/O3 are almost independent of the strain. Ac- 
cordingly, the longitudinal in-plane AMR is almost con- 
stant and the corresponding curves of piong are nearly 
identical (not shown). It should be noted again that for 
all samples under study the condition |p3| ^ |pi| does 
not apply and thus Eq. (|32p fails to correctly describe 
the angular dependence of p\ong- 

Whereas for configuration I (see above) and configu- 
ration III (not shown) the overall angular dependence 
of Piong is basically the same in the whole range of Szz 
under investigation, an essential variation with the strain 
occurs when H and thus M is rotated perpendicular to j 
(configuration II). In Figs. [TOTb) and (c) this variation is 
shown for j \\ [100] and j \\ [110], respectively. The com- 
plex angular dependences emerging for Ezz > —0.46%, 
i.e., the appearance of more than two maxima, qualita- 
tively agree with the polar plot depicted in Fig. [H They 
result from a competition between the second-order term 
P2Tn^ and the fourth-order term P4m^, and occur when- 
ever p2 and p4 are of comparable magnitude and opposite 
sign, or briefly, whenever P2+P4 is close to zero. As shown 
by the dashed curves in Figs. [Hl^a) and (b), this apphes 
for Ezz > —0.46% in agreement with Figs. fTOf b) and (c). 
For both J || [100] and j \\ [110],p2 + P4 monotonously in- 
creases with Ezz and changes sign at the transition from 
tensile to compressive strain. In the case of j \\ [100] 
symmetry demands that p2 + Pi exactly equals zero at 
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FIG. 10: (Color online) Variation of the angle-dependent nor- 
malized longitudinal resistivity {piong — po)/po with vertical 
strain e^^ at poH = 0.65 T and T = 4.2 K. The resistivity 
was recorded as a function of magnetic field orientation for 
(a) j II [110] and H in the (001) plane, (b) j \\ [100] and 
H in the (100) plane, and (c) j \\ [110] and H in the (110) 
plane. The experimental data are depicted by red circles and 
the calculated curves by black solid lines. 



Ezz = 0. The corresponding polar plot is schematically 
illustrated by the dash-dotted line in Fig. [5] 

The evolution of the transverse resistivity with Ezz is 
exemplarily discussed in Fig. [TT] where the angular de- 
pendence of ptrans/po is dcpictcd for j II [110] and m ± j 
(configuration II). For magnetization orientation along 
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FIG. 11: (Color online) Variation of the angle-dependent nor- 
malized transverse resistivity ptrans/po with vertical strain 
e^z at noH = 0.65 T and T = 4.2 K. The resistivity was 
recorded as a function of magnetic field orientation for j \\ 
[110] and H in the (110) plane. The dashed line was calcu- 
lated setting all magnetic anisotropy parameters to zero and 
demonstrates the influence of the MA at 0.65 T for Ezz = 
0.22%. 

[GDI], i.e., perpendicular to the layer plane, /Otrans equals 
P6 + Ps- As demonstrated in Fig. Efc), the sum pg + Ps 
is identical for both current directions and linearly de- 
creases with increasing Ezz- In Fig.jHJa) the dependences 
of = —{B + h) and = —{D-\-di) on e^z are dis- 
played as separate curves. They reveal that the values of 
P6 and ps have opposite sign and differ by more than one 



order of magnitude. The dominant parameter pg corre- 
sponds to the linear term m„ in Eq. (|19p which is usu- 
ally associated with the AHE. The increasing deviation 
of the curves in Fig. [TT] from a cosinusoidal oscillation 
for increasing compressive and tensile strain is primarily 
due to the influence of the magnetic anisotropy even at 
0.65 T. For comparison, the resistivity curve calculated 
for vanishing MA and £zz = 0.22% is drawn as a dashed 
line. 



V. SUIMIMARY 

Based on the expansion of the resistivity tensor with 
respect to the direction cosines of the magnetization up to 
the fourth order, we have presented a macroscopic analyt- 
ical model for the longitudinal and transverse resistivies 
of single-crystalline cubic and tetragonal ferromagnets. 
The model applies to arbitrary magnetization orienta- 
tions and current directions. It correctly describes the 
results of angle-dependent magnetotransport measure- 
ments performed on a series of (Ga,Mn)As layers with 
a vertical strain gradually varied from tensile to com- 
pressive. The resistivity parameters, obtained by fitting 
Eqs. (|18p and p9p to the experimental data, were found 
to systematically vary with the strain. 
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